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ABSTRACT 

Knowledge of the intrinsic shape of galaxy clusters is very important in investigating cos- 
mic structure formation and astrophysical processes. The reconstruction of the 3-dimensional 
structure usually relies on deprojecting 2-dimensional X-ray, Sunyaev-Zeldovich (SZ) and/or 
gravitational lensing observations. As known, a joint analysis of these data sets can provide 
the elongation of the cluster along the line of sight together with its length and width in the 
plane of the sky. An unbiased measurement of the Hubble constant can be also inferred. Due 
to some intrinsic degeneracies, the observational constraints obtained from such projected 
data-sets are not enough to allow an unique inversion. In general, the projected maps can be 
at the same time compatible with prolate, oblate and with many triaxial configurations. Even 
a prolate cluster might be interpreted as an oblate system and vice versa. Assuming that the 
cluster is axially symmetric is likely to overestimate the intrinsic ellipticity, whereas the sys- 
tem always looks rounder performing the inversion under the hypothesis of a triaxial cluster 
aligned with the line of sight. In general, analysing triaxial clusters under the prolate or oblate 
assumption may introduce strong biases even when the clusters are actually near to axial sym- 
metry whereas the systematics introduced assuming the cluster to be aligned with the line of 
sight are more under control. 
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- distance scale - gravitational lensing - cosmic microwave background 



1 INTRODUCTION 



The determination of the intrinsic shape of astronomical objects 
is a classic topic. The intrinsic structure of galaxies or cluster of 
galaxies directly probes the cosmic structure formation su ggest- 
ing h ow material aggregates from large-scale perturbations dWestl 
1 19941) . It also contains evidence about the nature and mechanisms 
of interaction of baryons and dark matter, since processes such as 
virialization, dissipation or gas cooling te nd to make systems mo re 
spherical, especially in the inner regions dKazantzidis et alj|2004h . 

The comp lex structure of halos also affects mass estimates 
dGavazzill2005l) and could cause a significant bias in estimating the 
inner den sity matter slope an d the concentration parameter of mat- 
ter halos dOguri et al. 2005). These quantities are crucial for any 
attempt at high precision cosmology and when comparing obs erva- 
tions with theoretical predictions from numerical simulations dVoitl 
120051) . 

The first attempts to determine three dimensional morpholo- 
gies were based on statistical approache s consist i ng in t he inversion 
of the distribution of apparent shapes. iHubblel (l926) first deter- 
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mined the relative frequencies with which galaxies of a given intrin- 
sic ellipticity, oriented at random, are observed as having various 
apparent projected ellipticities. Several following studies have then 
applie d similar methods to different classes of astronomical ob- 
jects dNoerdlingerjl 1979 : Bingeeli 198p|; Birmev & de Vaucquleursl 
198ll;lFasano & Vidll99ll:lde Theiie et alJll995l:lMohr et alj|l995l; 
Basilak os et al. 2000; Cooray 2000; Thakur & Chakrabortvil200ll; 



Alam & Ryder] |2002| ; iRvderil 1 19961 : IPlionis et al.l [2004 IPaz et al . 
2006). With the exception of disc galaxies, prolate-like shapes ap- 



pear to dominate all cosmic structure on a large scale. 

Together with statistical studies, the deprojection of single ob- 
jects wa s also investigated. Based on the Fourier slice theorem, 
Rybicki (1987) argued that, due to a cone of ignorance in the 
Fourier space, th e deprojection can not be unique, even assuming 
axial symmetry. iGerhard & Binnevl d 19961) further showed that ei- 
ther discy or boxy 3-D distributions can be compatible with the 
same projected image. 

Clusters of galaxies have the strong advantage that their struc- 
ture can be routinely probed with very heterogeneous data-sets 
at very different wave-lengths. This consideration boosted stud- 
ies on how combining X-ray surface brightness and spectral ob- 
servations of the intra-cluster medium (ICM), Sunyaev-Zeldovich 



2 M. Sereno 



effect (SZE) and gravitational lensing (GL) observations. One of 
the main result that can be obtained combining different data-sets 
is breaking the degeneracy between the line of sight elongation and 
the distance to the cluster, whic h allows in princi ple an unbiased 
estimate of the Hubble constant dFox & Penll2003) . The deprojec- 
tionof the density distribution was also reconsidered. IZaroubi et al.l 
d 19981) proposed a method based on the axial symmetry assumption 
and on the extrapolation of the image Fourier transform into the 
cone of ignorance. Alternative procedures can be based either on 
the iterative Ric hardson Lucy de c onvolution and again assuming 
axial symmetry jReblinskvlExx]: Puchwein & Bartelmann 2006) 
or on a per turbation approach ([Pore et alJl200lh . lDe Filippis et al.l 
d20Q5h and lSerenoetaU {2006) finally applied a joint X-ray plus 
SZE parametric analysis to a sample of 25 clusters, finding that 
prolate rather than oblate shapes seem to be preferred, with signs 
of a more general triaxial morphology. 

Many of the above methods, despite very insightful, often rely 
on very restrictive assumptions, sometimes not clearly stated, and 
aim to obtain very general results. The methods are then usually 
tested with the application to some numerical simulations, so that 
possible degeneracies in the deprojection techniques can be eas- 
ily over-sought. In this paper I propose a simple analytical discus- 
sion aimed to explore what are the features of the cluster shape 
that we can really measure by combining X-ray, SZE and gravi- 
tational lensing data and under which conditions the deprojection 
can be performed unequivocally. The paper is as follows. In Sec. [2] 
I describe how a 3-D ellipsoid casts on the plane of the sky. Sec- 
tion [3] briefly presents the main characteristics of X-ray, SZE and 
GL observations, whereas Sec.|4]discusses which constraints about 
cluster shape and orientation can be inferred. In Sec. [5] I consider 
when only a single deprojection is allowed, whereas Sec. [6] treats 
the error made on the determination of the cluster shape under par- 
ticular assumptions. Section|7]is devoted to when and how a prolate 
cluster can appear as oblate and vice-versa. At the end, Section [8] 
contains some final considerations. 



2 PROJECTED ELLIPSOIDS 

High resolution TV-body simulations have clearly shown that the 
density profiles of matter halos are aspherical and how such pro- 
files can be accur ately described by concentric triaxial ellipsoids 
with aligned axes Uing & Sutol2002h . Then, the electron density of 
the intra-cluster medium can be assumed to be constant on a family 
of similar, concentric, coaxial ellipsoids. The intra-cluster medium 
(ICM) distribution in clusters of galaxies in hydrostatic equilibrium 
traces the gravitational potential. Since we are considering a triaxial 
elliptical gas distribution, the gravitational potential turns out to be 
constant on a family of similar, concentric, coaxial ellipsoids. El- 
liptical gravitational potential can turn unphysical for extreme axial 
ratios, giving negative density regions or very unlikely configura- 
tions, but as far as inner regions are considered, they can provide 
very suitable approximations. 

In an intrinsic orthogonal coordinate system centred on the 
cluster's barycentre and whose coordinates are aligned with its 
principal axes, a spheroidal ICM profile can be described by only 
one radial variable £, 



t 2 = 2 2 



(1) 



in units of a scale-length r c . Without loss of generality, we can fix 
ez = 1, which means that our reference scale-length is the core ra- 
dius along S3. i n t. For an axially symmetric cluster, if the polar axis 
is aligned with the third coordinate axis, ej = {1/f/mt, l/^int, 1} 
for a prolate model and = {gint, Qint, 1} for an oblate model, 
where gi n t ^ 1 is the intrinsic axial ratio. 

When viewed from an arbitrary direction, quantities co nstant 
on sim ilar ellipsoids project themselves on similar ellipses dStarkl 
[l977h . Three rotation angles relate the intrinsic to the observer's 
coordinate system, i.e. the three Euler's angles, 9eu, fEu and ipEu 
of the three principal cluster axes with respect to the observer. A 
rotation through the first two Euler's angles is sufficient to align 
the i3 i0 b s -axis of the observer coordinates system {a^obs} with 
whatever direction. In what follows, I will assume that the a^obs- 
axis is aligned with the line of sight to the observer, i.e. the direction 
connecting the observer to the cluster centre. Then, in the intrinsic 
system, the line of sight has polar angles {8, <j>} = {9e u , fEu — 
7r/2}. With a third rotation, ^eu, we can properly align the a;i, bs- 
and a;2,obs-axis in the plane of the sky. In general, ipEu is the angle 
in the plane of the sky between the projection of the a^mt-axis and 
the £2,ob s -axis. Unfortunately, the direction of the a^jnt-axis is not 
known, so that we can not choose a reference system in the plane 
of the sky such that tpEu = 0. Then, if not stated otherwise we will 
line up the a7i j0 bs- and £2,obs-axis with the axes of the projected 
ellipses. 

The ellipticity and the orientation of the projected ellipses de- 
pend only on the intrinsic geometry and orientation of the system. 
The axial ratio of the major to the minor axi s of the observe d pro- 
jected isophotes, e p (^ 1), can be written as jBinggelill 19801) . 



3 + 1+ V(j-0 2 + 4fc 2 

6p \ j + i- Vti - 2 + 4fc 2 ' 

where j, k and / are defined as 

J 



2 2 2 2 2 2 

= e 1 e 2 sin #e u + e x cos feu cos <^e u 



+ 62 cos 2 6eu sin 2 </pe u , 
k — (e 2 — e 2 ) sin^EuCosipEuCosfeu, 
I = e 2 sin 2 v?Eu + e2 cos 2 v?e u . 



(2) 



(3) 
(4) 
(5) 



As written before, the projected direction of the a?3,i nt -axis is 
not known. If we assume that the coordinate a xes in the pla ne of 
the sky lie along the axes of the isophotes then l lBinnevll 19851) 



ipEu — g ar ctan 



2k 



(6) 



The apparent principal axis that lies furthest from the projection of 
th e Z3,int-axis o nto the plane of the sky is the apparent major axis 
if teinnevll 19851) 



(j — I) cos 2-ipEu + 2k sin 2i/>eu ^ 



(7) 



or the apparent minor axis otherwise. In general the right-hand side 
of Eq. (JSJl is the angle between the principal axes of the observed 
ellipses and the projection onto the sky of the a^mt-axis. 

The observed cluster angular core radius 9 P is the projection 
on the plane of the sky of the cluster angular intrinsic core radius 
dStarkll 19771) , 



(8) 



Along each axis, et is the inverse of the corresponding core radius 



where 8 C = r c /Dd, with the angular diameter distance to the 
cluster, and / is a function of the cluster shape and orientation, 
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' = e\ sin 2 6>bu sin 2 y> E u + e 2 sin 2 # Eu cos 2 <^e u + cos 2 Eu . (9) 



3 MULTI-WAVELENGTH OBSERVATIONS 

In this section, I briefly summarise the main features of X-ray, SZE 
and GL observations. As a reference model, we consider a cluster 
electron density distribution, n e , described by an ellipsoidal triaxial 
/3-model. In the intrinsic coordinate system, 



n c = n e o [ 1 + — 



-3/3/2 



(10) 



where /3 is the slope, r c is the core radius and n c o is the central 
electron density. We are interested in observed quantities which are 
given by projection along the line of sight of powers of the electron 
density n c . The projection of a generic power m o f n e , as given in 
Eq. {TO}, 

can be written as l lDe Filippis etailfe oOS ) 



l (l)dl 



n c0 V7r- 



-T [3m/3/2 - 1/2] D d 6 c 
r [3mP/2] ^/ 

(l-3mfl/2 



=p"2 



(ID 



'P 

1 

eA 



(12) 



where 9i = Xi t0 h s /Dd is the projected angular position on the 
plane of the sky of Xi.obs- The quantity r c /y/J in Eq. dl It repre- 
sents the half-size of the ellipsoid along the line of sight. It can be 
conveniently rewritten as 

Vf 

which represents the definition of the elongation eA; r p (= D d 9 p ) 
is the projected radius.If ca < 1, then the cluster is more elongated 
along the line of sight than wide in the plane of the sky. In terms of 
the elongation, the projected core radius can be expressed as 

e^eJ^Y' 3 . d3) 

V eie2 / 

In a Friedmann-Lemaitre-Robertson-Walker universe filled 
with pressure-less matter and with a cosmological constant, the an- 
gu lar diameter distan ce between the redshift z d and a source at z B 



is dSereno et al ]|200ll . and references therein) 



£>0d, 



with 



£(z) 



1 
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Ho 1 + z s If2p 



-Sinn 



£(z 



-dz 



H( 



Ho 



(14) 



(15) 



j(l + z) 3 + Qao + Oko(1 + z) 



where Ho, Qmo and S7ao are the Hubble parameter, the normalised 
energy density of pressure-less matter and the reduced cosmolog- 
ical constant at z = 0, respectively. 57ko is given by £Iko = 
1 — CImo — Oao, and Sinn is defined as being sinh when Qko > 0, 
sin when fijco < 0, and as the identity when Q,ko = 0. 



3.1 X-ray Surface Brightness 

Cluster X-ray emission is due to bremsstrahlung and line radia- 
tion resulting from electron-ion collisions in the high temperature 
plasma (ksTc ~ 8-10 keV, with fce being the Boltzmann con- 
stant). The X-ray surface brightness Sx can be written as 



Sx 



4tt(1 + z) 4 



n c A E (T c ,Z)dl, 



(16) 



where A c is the cooling function of the ICM in the cluster rest frame 
and depends on the ICM temperature T c and metallicity Z. Assum- 
ing an isothermal plasma with constant metallicity and taking the 
result from Eq. Jilt for m = 2, we get 



Sx = Sxo 1 



q2 , c 2fl2\ 1/2-3/3 

h + e P y 2 



where the central surface brightness Sxo reads 

A c 



Sxo 



2 D d e p r(3/3 - 1/2) 

f n eO" 



40F(l + z)4 cu e A r(3/3) 



(17) 



(18) 



3.2 The Sunyaev-Zeldovich Effect 

Photons of the cosmic microwave background (CMB) that pass 
through the hot ICM of a cluster interact with its energetic 
electrons through inverse Compton scattering, slightly distort- 
ing the CMB spectrum. This is the Sunyaev-Zeldo vich effect 
(SZE) JSunvaev & ZeldovicrJ [l97Cl: iBirkinshawl 1 19991) . which is 
proportional to the electron pressure integrated along the line of 
sight. The measured temperature decrement ATsz of the CMB is 
given by, 



Tcm: 



= fsz{v,T e 



, OT^E 



n e T e dl, 



(19) 



where Tcmb is the temperature of the CMB, <tt the Thompson 
cross section, m c the electron mass, c the speed of light in vacuum 
and fsz(v, T e ) accounts for relativistic corrections at frequency v. 
For an isothermal /3-model, taking the result in Eq. i ll It for m = 1, 
we obtain, 



ATsz = AT 1 + 



?i + Cp02 2 
01 



1/2-3/3/2 



(20) 



where ATo is the central temperature decrement which includes all 
the physical constants and the terms resulting from the line of sight 
integration 



ATo 



-LCMBJSZ(V, J c ) — n c o\ 



D d 9 p T(3/3/2 



1/2) 



EA 



T(3/3/2) 



(21) 



3.3 Gravitational Lensing 

Clusters of galaxies act as lenses deflecting light rays from back- 
ground galaxies. In contrast to SZE and X-ray emission, gravita- 
tional lensing does not probe the ICM but maps the total mass. The 
cluster total mass can be related to its gas distribution if the intra- 
cluster gas is assumed to be in hydrostatic equilibrium in the cluster 
gravitational potential. If we assume that the gas is isothermal and 
that non-thermal processes do not to contribute significantly to the 
gas pressure, the total dynamical mass density can be expressed as 



ptot 



.-(. 



V 2 (lnn c ), 



(22) 



where G is the gravitational constant and /im p is the mean particle 
mass of the gas. Ellipsoidal ICM distributions determine ellipsoidal 

otentials, which are w idely used in gravitational lensing analyses 

Schneider etalj 1992b . 
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If we assume that the ICM follows a /3-model distribu- 
tion, the projected ma ss density can be subsequently derived 
foe Filippis et al.ll2005h . The lensing effect is determined by the 
convergence k — E/E cr which is the cluster surface mass density 
in units of the critical density E cr , 

c 2 D s 



— " AtvG D d D ds ' (23) 

where D dB is the angular diameter distance from the lens to the 
source and D s is the angular diam eter distances from the o bserver 
to the lens. The convergence reads (De Fili ppis et al . 2005) 

-3/2 



k = k 1 



where 



?! + el 



l+e p 



1 + 



9i + egflj 

el 



3n/3- 



fc B T c 



1 n , 2s 1 D is 



(24) 



(25) 



r/imp eA 

The convergence can be measured through either a weak lensing 
analysis of statistical distortion of images of background galax- 
ies or by fitting the observed surface mass density to multiple im- 
age strong lensing systems. Although the hypotheses of hydrostatic 
equilibrium and isothermal gas are very strong, total mass densities 
obtained under such assumptions can yield accurate estimates even 
in dynamically active c lusters with irregular X-ray morphologies 
dDe Filippis et alj|2004l) . 



4 JOINT ANALYSIS OF HETEROGENEOUS DATA SETS 

This section discusses which geometrical constraints on the cluster 
shape can be inferred from projected maps. 



4.1 Unknown parameters 

The intrinsic shape of an ellipsoidal cluster is described by two axis 
ratios, ei and e2. Its orientation is fixed by three Euler's angles, 
0e u , </>Eu and ipEu- 

The density ICM profile is characterised by some other scale 
parameters. For a /3-model, three additional parameters come in: 
the central density normalisation, n c o, the slope index, and the 
core radius r c . Then, 5 + 3 parameters characterise the ICM distri- 
bution. Under the hypothesis of isothermality, a single value, T e , 
characterises the temperature of the cluster. If the metallicity is 
nearly constant, then one other parameter Z is enough to describe 
this quantity. 

The cosmological dependence enters through the cosmologi- 
cal distances. For a flat model of universe, the distance-redshift re- 
lation is determined by two parameters: the Hubble constant, Ho, 
and the matter density parameter Qmo- 



4.2 Observational constraints 

As seen before, 3-D ellipsoids are casted in 2-D projected ellipses. 
By fitting an elliptical profile to the X-ray and/or SZE data, both the 
projected axis ratio, e p , and the orientation angle can be measured 
providing two constraints on the intrinsic shape 

e p = e p (ei,e 2 ;6 l Eu,i?Eu), (26) 
ifeu = </>Eu(ei,e 2 ;6>Eu, i)eu). (27) 

Eq. j27\ expresses the freedom to align the isophotes with the co- 
ordinate axis. 



Together with the shape and orientation of the ellipses, the 
fitting procedure can also provide further constraints on the density 
profile. In particular, for a /3-model, the slope /? and the projected 
core radius P can also be determined from data. 

/3obs = P, (28) 
P = 6 l p (ei,e 2 ;eA(ei,e2;6'Eu,i?Eu)), (29) 

where we have made clear the dependence on the unknown quanti- 
ties. 

Besides these constraints on the density profile, the temper- 
ature of the ICM, as well as its metallicity, can be inferred from 
spectroscopic X-ray observations with sufficient spectral resolution 



T c = 
Z = 



2>obs ■ 



(30) 
(31) 



The observed values of the central surface brightness, Sxo, 
Eq. d!8t , of the central temperature decrement, ATq, Eq. d21t . and 
of the gravitational lensing convergence, Eq. ( 125b . provide three 
further constraints, 



2 

Sxo oc — —D d , 

eA 

Al oc D d , 

eA 

, 1 D Ab 
k oc — . 

eA -Ds 



(32) 
(33) 
(34) 



The dependence on the elongation of the cluster shows up both 
in the relation between the projected and intrinsic core radius, 
Eq. l |29t , and in the expression of the central quantities, Eqs. 1 1321 
1341 . The convergence k depends on the cosmology through the ratio 
of distances Dd B /D s . Therefore k depends only the cosmological 
density parameters Qj, and not on the Hubble constant Ho. In all, 
we have two observational constraints less than the unknown pa- 
rameters, so that the system is under-constrained. 

4.3 Inferred quantities 

Let us see what we can learn on the cluster structure. Equation l !29t 
and Eqs. J30II3U refers specifically to an isothermal /3-model with 
constant metallicity. In any case, they show a general feature: the 
parameters which describe the temperature and density profile can 
be derived in principle with accurate spectroscopic and photomet- 
ric observations even for more complicated models. Instead, we are 
mainly interested in the intrinsic shape (2 parameters) and orienta- 
tion (3 parameters) of the cluster. As a first step, let us assume that 
the cosmological density parameters are independently known, so 
that cosmological distances are known apart from an overall factor 
proportional to the Hubble constant, Di oc c/ Ho - Then, the system 
of Eqs. d32|34fc is closed and can be simply solved, 

Sxo 
ATb' 

' (36) 



n c o 



(35) 



f'A 



Ho 



ko 
k 



Sxo 
A7f 



(37) 



As well known, combining SZ and X-ray observations only fixes 
the central gas density. On the other hand, the degeneracy between 
the Hubble constant and the physical size of the cluster along the 
line of sight can be broken only with the additional information pro- 
vided by gra vitational lensin g under the assumption of hydrostatic 



vided by gravitational lensing t 
equilibrium iFox & Penll2002f) . 
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Observations of either multiple strong lensing image systems 
or weak lensing of background sources with well determined pho- 
tometric redshifts could provide further information on the cosmo- 
logical parameters fij, In fact, the value of fco changes according 
to the redshift z s of the lensed source through the ratio Dds/Ds, so 
that for each background source redshift we have an additional con- 
straint on the geometry of the universe. With a sufficient number of 
image systems at different redshifts an estimate of al l cosmologica l 
parameters involved could therefore be performed (Sereno 2002). 
Stacking gravitational lensing data of systems at different redshifts 
will then provide information on the cosmology and will reduce 
the observational uncertainties on the density parameters but will 
not add constraints on the intrinsic shape of the cluster. One im- 
age system is enough to break the degeneracy between the Hubble 
constant and the elongation of the cluster. On the other hand, if we 
trust independent estimates of the Hubble constant, then SZ and X- 
ray data are enough to measure the elongation of the cluster and 
we do not need the additional theoretical constraint of hydrostatic 
equilibrium and the observational gravitational lensing data. 

From the above analysis it is clear that even relying on extraor- 
dinarily accurate projected X-ray, SZE and GL maps, the intrinsic 
shape of the cluster can not be unambiguously inferred. The only 
information that we can establish on the cluster is its width and 
length in the plane of the sky, i.e. r p and r p /e p , and its size along 
the line of sight, r p /eA; r p is related to the intrinsic scale-length 
through Eq. l |29l l. Furthermore, we have a relation which expresses 
the third Euler's angle, ip^u in terms of the other intrinsic parame- 
ters. These 4 constraints are what an astronomer dealing with pro- 
jected X-ray, SZE and GL maps can use to determine the (2) in- 
trinsic axial ratios, the (3) orientation angles and the (1) intrinsic 
length-scale of the cluster. 

4.4 On a general density distribution 

The i sothermal j3 model can be sometimes inaccurate iRasia et all 
2006) but the conclusions on what we can learn on the 3-D struc- 
ture of galaxy clusters from projected maps do not change con- 
sidering more accurate models. The temperature profile should ac- 
count for a central cool re gion, if any, and a gradient al large 
radii ( Vikhlini n et al. 2006). In the same way, the ICM density 
profile should be able to account for a possible central steep in- 
crease of the surface b rightness and a change of slope at larger radii 
dVikhlinin et al.l2006h . Accurate modelling is required when study- 
ing the ICM physics or when reliable mass estimates have to be 
obtained from deep observations but, as far as an insight on which 
physical constraints on the 3-D shape we can get from projected 
maps is concerned, the conclusions are not affected. Constraints 
obtained using more accurate models are more reliable but do not 
allow to break any degeneracy discussed in the framework of the 
(3 model. In fact degeneracies are connected to the ellipsoidal in- 
trinsic structure, not to the specific features of the radial profile. 
Then, whatever the model used, what can be inferred on the clus- 
ter intrinsic structure is its elongation, ellipticity and orientation of 
the isophotes, together with an estimate of the Hubble constant (if 
GL data are available) and with the parameters characterising the 
properties of the density distribution (slope, concentration,...). 

Under very broad assumptions, i.e. ICM density profile mono- 
tonically decreasing with radius, any radial distribution (tempera- 
ture, metallicity,...) can be expressed as a function of the ICM dis- 
tribution (or of the gravitational potential). Furthermore, with the 
hypothesis of hydrostatic equilibrium, the temperature profile is de- 
scribed by the same axis ratios and orientation of the ICM. Then, 



both the X-ray surface brightness, Eq. J17t . and the SZ tempera- 
ture decrement, Eq. (19), are projections of some function of n e . 
Furthermore, the temperature profile can be in principle extracted 
taking ratios of the X-ray emission in differ ent energy band s and is 
therefore independent of the elongation dFox & Pen 2002). Then, 
a more detailed temperature analysis would constrain important 
features of the temperature profile, such as a possible cool core 
radius, but it would not help in providing additional independent 
constraints on the shape and orientation of the cluster. 

In general, an intrinsic volume densit y F v and it s projection 
on the plane of the sky F s are related to by dStarklll977h 

•"o-Tfjf^TF^* <38) 

where £ is the ellipsoidal radius in the plane of the sky. For coordi- 
nate axes oriented along the isophotes, £ 2 = D\(0\ + ep^O 2 ) f / 
£ 2 can be rewritten in terms of generic intrinsic (l c ) and projected 
(l p ) scale-lengths, related as in Eq. C2J, i.e. £ 2 = D\(0l + 
e p@2)(lc/l P ) 2 ■ Since the integral in £ in Eq. d38b is proportional 
to the intrinsic scale length, we can write 

Q 

F s oc Da— f a (ep,ip Eu ;0 p ;pi,... ) (39) 
eA 

where as usual l p = Dd0 p and pi are the parameters describing 
the 3-D density function F v , The structure of Eq. l |39t is the same 
of Eqs. < |32b and l |33t , which were obtained for a specific model. 
The dependence on the elongation ba is decoupled from the depen- 
dence on the apparent ellipticity and inclination and the parameters 
characterising the 3-D profile only accounts for the radial depen- 
dence of the projected density. Then, as f3 can be derived by fitting 
the projected isophotes to a /3-model, the parameters of a differ- 
ent profile (slope, concentration, ...) can be determined by a similar 
procedure as well. As before, the intrinsic scale-lengths are related 
to their projected values as in Eq. d 12b . So the statement that the 
only information that we can establish on the cluster is its sizes in 
the plane of the sky and along the line of sight does not depend on 
the particular mass model of the cluster. 

The above considerations show that luminosity and surface 
brightness observations in the optical band (which are other exam- 
ple of projected maps) would provide the same kind of information 
on the intrinsic structure as those obtained by X-ray and SZE ob- 
servations. 



5 PARTICULAR SOLUTIONS 

Having ascertained that a full inversion can not be performed only 
based on 2-D maps, I now turn on the possibility of deproject- 
ing cluster observations under particular assumptions. Let us first 
consider some intrinsic degeneracies. As well known, when deal- 
ing with projected maps, we are not able to determine which ex- 
tremity of the cluster is pointing towards the observer, e;(#Eu) = 
ei(— #Eu), i = {p, A}. That is why in what follows we will limit 
our considerations to ^ 6>e u ^ n/2. In the same way, the ob- 
served quantities are invariant for i^eu — > <£e u + 7T. Elongation 
and projected ellipticity are also symmetric around ip^u ~ 7r/2, 
i.e. eA(7r/2 + ^ Eu ) = 6a(7t/2 — <£e u ) and the same for e p . Fur- 
thermore for ip-Eu = tt/4, e p and eA do not change for ei — > e-z 
and e2 — ► ei. 

One simple way to break intrinsic degeneracies is to fix the 
orientation of the line of sight in the intrinsic system of the ellip- 
soid. This lowers the number of unknown parameters by bringing 
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Figure 1. Solution domains as a function of the orientation angles of the line of sight for three different triaxial ellipsoids. The grey, black and white regions 
denote the loci of line of sight directions in which only the prolate, only the oblate and both solutions are admissible, respectively, is fixed to 1. Left panel: 
a nearly prolate ellipsoid, ei = 1.3, e 2 = 1.4. Middle panel: a nearly oblate ellipsoid, e\ = 0.7, e 2 = 0.6. Right panel: a pretty triaxial ellipsoid, ei = 1.2, 
e 2 = 0.8. 
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Figure 2. Solution domains as a function of the intrinsic axis ratio. e$ is 
fixed to 1. The white and black spotted regions denote the loci of axial ratio 
for which the double and the oblate solution are more likely, respectively. 



out c9eu and </?eu- Since observations can determine the elongation 
and e p and since the third Euler's angle can be expressed as func- 
tion of the axial ratios and the other two angles, one could think that 
assuming that #e u and c^eu are known, then the inversion to deter- 
mine e\ and e 2 could be unambiguously performed. Unfortunately, 
the equations for e\ and e 2 



eA(ei, e 2 ; 



VEu| flx ) = e A | ob 



e P (ei,e 2 ; 0Eu| flx , VEu| fix ) = e P 



(40) 
(41) 



are not linear. Multiple solutions can exist even when the orienta- 
tion of the system is correctly known. The situation gets even worse 
when the orientation is not exactly fixed to the actual direction. 
When performing a deprojection, one can make an attempt to fix the 



inclination angles to some particular (trial) values and then check 
if the inversion is possible under this assumption. Unfortunately, 
choosing a trial orientation for the line of sight we are not even as- 
sured about the existence of even one solution. Just as an example, 
if the line of sight of a triaxial system with ei = 0.7 and e 2 = 1.2 
lies along cos feu = 1/2 and c/peu = tt/3, then the equation sys- 
tem for ei and e 2 does not have a real solution even if the trial val- 
ues are fixed to cos0E U | trial = 1/2 and ( ^Eu| trial — ¥>Eu)/¥>Eu = 
10 -2 , i.e with a very small error Ac/?eu ~ 0.5 deg. On the other 
hand, for (v?Eu| trial — veiO/veu = — 10~ 2 , there are two solu- 
tions. 

It can then be useful to consider some particular assumptions 
under which the solution exists and is unique. Interesting config- 
urations are either an aligned triaxial or an axially symmetric el- 
lipsoid. Under the assumption that one of the principal axis of the 
ellipsoid is aligned with the line of sight, the inversion can be easily 
performed and the intrinsic axial ratios are easily r ecovered. This 
hypothesis was exploited in De Filipp is et alj d2005t) to deproject a 
sample of luminous X-ray clusters with SZE observations. As an 
example, when the a:3,i nt -axis is aligned with the line of sight, i.e. 
8eu ~ and i^eu = n/2, then the maximum between ei and e 2 
will be e p /eA and the minimum 1/eA (as usual e$ — 1). 

Another very popular ch oice is assuming tha t the cluster shape 
is nearly axially symmetric dSereno et al. 2006). In this case the 
ICM distribution is characterised by just two parameters: the ratio 
of the minor to the major axis, qi nt 1), and the inclination angle, 
i, between the line of sight and the polar axis. If the polar axis 
lies along the X3,int-axis, then i = #e u . The major (minor) axis 
of the isophotes coincides with the projection of the polar angle if 
the ellipsoid is prolate (oblate). A useful parameter to quantify the 
triaxiality degree of an ellipsoid is 



T = 



^mid 



(42) 



where e m t n , e m id and e max are ei, e 2 and sorted in growing 
order. Oblate and prolate clusters correspond to T = and 1, re- 
spectively. 

The projected ellipticity and the elongation of a prolate clus- 
ter can be easily expressed in terms of the intrinsic parameters 
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dSereno et ai1l2006l) . 



eA 



yi-(i-gL)cos 2 » 

9int 

1 - (l - gL) COS 2 » 
9int 



(43) 



(44) 



The previous couple of equations can be then easily inverted to 
infer the intrinsic shape, 



<7int 



cos i = e t 



(45) 



(46) 



A prolate-like solution is then admissible only when the size along 
the line of sight is larger then the minimum width in the plane of 
the sky, i.e. when 



ga ^ e p . 



(47) 



The relations between the intrin sic parameters of an oblate 
cluster and its observable features are dSereno et all2006T) 



(48) 



V i 1 ~ lint) COs2 i + ll 



eA = 
Then, 

<7int = 

cos i — 




— 1 cos 2 i. 



e p eA 



(49) 



(50) 



(51) 



An oblate-like solution is admissible only when the size along the 
line of sight is larger than the maximum size in the plane of the sky, 
i.e. when 



eA > 1. 



(52) 



Both the prolate and the oblate solutions are admissible at the same 
time only when 



1 ^ eA ^ e p , 



(53) 



i.e. when the size along the line of sight is intermediate with respect 
to the projected dimensions. Even the assumption of axial symme- 
try is not enough in general to have an unique solution, but one have 
to specify if the cluster is assumed to be either prolate or oblate. 

Let us now consider if, given a random orientation of the line 
of sight, it is more likely that the observed elongation and elliptic- 
ity are compatible with either a prolate or an oblate solution or with 
both of them. If clusters are randomly oriented, then they uniformly 
occupy the cos #eu - V?e u plane. A cluster will be more likely in- 
terpreted as prolate than as oblate if the total area of the loci in 
the cos #eu-<£eu plane where only the prolate solution is admissi- 
ble (grey regions in Fig.QJ is larger than that corresponding to the 
oblate case (black regions). In Fig. [TJ we show the regions in the 
cos feu-fEu plane where a prolate and/or oblate solution is pos- 
sible for three different sets of intrinsic axial ratio. You can also 
note the recurrence properties in the cos ^eu-^Eu- We consider a 
nearly prolate cluster with ei = 1.3, e2 = 1.4 (T = 0.75), a 
nearly oblate ellipsoid with ei = 0.7, e?. — 0.6 (T = 0.25) and a 
pure triaxial specimen, ei = 1.2, e2 = 0.8 (T = 0.5). Even when 



the cluster is intrinsically close to either a prolate (T = 0.75) or 
oblate (T = 0.25) geometry, projection effects nearly completely 
hide this property and conflicting domains in the cos#Eu - V^Eu 
plane have very similar extensions. However, for most of the orien- 
tations, both solutions are compatible with the observed quantities. 
These relative proportions change when the cluster get more triax- 
ial (T = 0.5). In that case the areas of the three different domains 
are similar. 

In Fig. [2] we consider the most likely domain as a function of 
the intrinsic axial ratios. The most likely condition, for fixed values 
of ei and e2, is determined by checking which existence conditions 
is fulfilled by elongation and projected ellipticity when averaged 
over the cos6» Eu - ^Eu plane, i.e. (eA) VEu ,e Eu > ( e p)<t> E u,e E u for 
a most likely oblate solution, (eA)v3 Eu ,e Eu < 1 for a most likely 
prolate solution and 1 < (eA)^ Eu ,e Eu ^ (e P }^ Eu ,e Eu for a most 
likely double solution. For a large range of axial ratios, the cluster 
projections are compatible with both solutions. This holds in par- 
ticular for the large central squared area in Fig. [2] corresponding 
to 0.4 < ei,e2 < 2.5, and in the region corresponding to nearly 
prolate clusters (ei ~ > 1). We see that in the parameter range 
explored in Fig. [2] 0.1 ^ ei, 62 ^ 10, there are no locations where 
the case of only prolate solution is the most likely. 



6 APPROXIMATE DEPROJECTIONS 

Let us now consider the error made when the intrinsic shape of the 
cluster if inferred under the assumption that the cluster is either 
aligned with the line of sight or is axially symmetric. Since we are 
not interested in observational noise but just in the intrinsic error 
due to degeneracies in the deprojection, we suppose to have fidu- 
cial X-ray, SZE and GL data and that we can analyse them using the 
proper density distributions. This implies that we can get a correct 
measurements of both elongation and ellipticity from observations. 
Elongation and ellipticity depend only on the intrinsic ellipticity 
parameters, ei, e2 and e^, and on the orientation of the line of sight, 
y?Eu and #e u , see Sec. [2] As discussed in Sec. [5] eA and e p are all 
we need to perform the deprojection under the hypotheses of either 
a triaxial system aligned with the line of sight or an ellipsoid of 
revolution. Then, in order to evaluate the error made in the depro- 
jection we have just to compare the intrinsic values (that we can fix 
at the beginning) to the ones inferred under the procedures outlined 
on Sec. [5] To this aim, it is useful to consider the maximum axial 
ratio, g max = e m i n /e max , i.e. the ratio of the minor to the major 
axis of the intrinsic ellipsoid. 

In Figs. [3] and [4] (upper panels), we consider the error made as 
a function of the line of sight orientation for two different sets of 
axial ratios: a nearly oblate case with e\ — 0.6, e2 = 0.7 (T = 
0.25, g m ax = 0.6) in Fig [3] and a triaxial case with ei = 1.2, 
e2 = 0.8 (T = 0.5, g max = 2/3) in Fig[4] The inversion under the 
respective hypothesis can not be performed for the directions filling 
the black regions. The situation for a nearly prolate cluster with 
T = 0.75 are very similar to the nearly oblate case with T = 0.25 
if we interchange the positions of the middle and of the right panel 
in Fig. [3] Under the hypothesis that the line of sight is along one 
of the principal axes of the ellipsoid, the axial ratio will be over- 
estimated, i.e. the cluster will end up looking rounder. Obviously, 
the error is minimum when the lining up is pretty much satisfied 
(9eu ~ 7r/2 and ^eu ~ or 7r/2, (9e u ~ 0). On the other hand, 
the loci of orientations where the deviation is maximum depend on 
the intrinsic geometry of the ellipsoid. For axial ratios close to the 
oblate case, see Fig[3] or close to the prolate geometry (T = 0.75), 
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Figure 3. Deviation of the inferred maximum axial ratio from the actual one (Aq m ax). The intrinsic parameters are fixed to e\ = 0.6, e2 = 0.7 and = f . 
<?max has been derived under the hypothesis of a triaxial ellipsoid aligned along the line of sight, of a prolate and of an oblate ellipsoid in the left, central 
and right panels, respectively. In the upper panels, Aq max is plotted as a function of the Euler's angle of the line of sight. Contours of equal derived q max 
are drawn in steps of 0.05. The thicks line tracess the loci of points where the deprojection gives the actual g max . The hypotheses breaks down in the black 
regions. In the bottom panels, the normalised probability distribution of the deviation Ag max is considered. Data are binned at intervals of 0.01. The shadowed 
and light shadowed regions contain the 68.3% and the whole range of the inferred values, respectively. 



we have Ag max > 0.25 for jpeu ~ tt/2 and cos feu ~ 0.6. For 
a T — 0.5 shape, see Fig[4] the deviation is maximum for ifEu ~ 
7r/4 and Eu < 7l"/2, where A<j max > 0.30. 

Under the hypothesis of axial symmetry, the maximum axial 
ratio can be either over- or under-estimated. Let us first consider 
the nearly oblate intrinsic shape represented in Fig. [3] Under the as- 
sumption of oblate geometry (right panel), for most of the orienta- 
tions the error will be in the range —0.1 < Aq max < 0.1; the max- 
imum deviation, Ag max < 0.25, occurs near the forbidden region 
for large values of </?e u (~ tt/2). Deviations will be much larger 
under the assumption of a prolate geometry (middle panel). It can 
be easily seen that the region where —0.1 < Ag max < 0.1 is much 
smaller than the corresponding one obtained under the oblateness 
assumption. For a triaxial intrinsic shape, Fig. [4] the oblate and the 
prolate hypotheses works pretty much in the same way. 

The normalised probability distributions of the deviation 
Ag max are plotted in the bottom panels of Figs. [3] and [4] The 
distributions have been obtained under the hypothesis of ran- 
domly oriented clusters, so that the probability of a deviation, 
P(Ag max )dAg max is proportional to the area in the <^e u -cos #e u 
plane of the region where the values of the deviation are between 
Ag max and Aq max + dA<j max . The portions of the grey regions 
above or below the averaged value refer to the same total area in 
the cos Oeu-^Peu plane. Both the light grey regions below and above 
the average correspond to a region in the cos #e u -v?e u plane whose 
area is 15.8% of the total. The fact that the light grey region be- 
low the average has a smaller extension than that above means that 



in the corresponding region in the cos #e u -(^eu plane the range of 
Ag max values is much smaller than in the other one. As general 
feaures, we see that the aligned triaxial hypothesis over-estimate 
the axis ratio and that, in general, there is always a long tail for 
large positive values of Ag max , i.e. there is always a chance that 
the cluster appears nearly round. If the case of T — 0.25, the av- 
erage value of the deviation is 0.09, —0.06 and < for an aligned 
triaxial, prolate or oblate assumption, respectively. 

The total range of inferred g max under different hypotheses 
and for various orientations of the line of sight depends in gen- 
eral on the triaxiality degree of the cluster. In Figs. [5] and [6] we 
show the general features of the inferred axial ratio as a function 
of the triaxiality parameter T for two different intrinsic axial ra- 
tios, </ max = 1/3 and 2/3, respectively. In particular, we consider 
the value of the inferred g max averaged over all the possible ori- 
entations of the line of sight, the range including the 68.3% of the 
inferred axial ratios around the averaged value (grey regions), and 
the total range of inferred values (light shadowed region). These 
values try to summerize the full information contained in the 9e u - 
ifEu plane, as shown for single T values in Figs.[3]and|4] 

The extension of the grey regions in Figs. [5] and [6] has the 
same meaning as in the bottom panels of Figs. [3] and [4] If we 
make the hypothesis of a triaxial cluster aligned with the line 
of sight, the error is nearly independent of the effective triaxi- 
ality. The axial ratio g max will be overestimated in average by 
0.07(0.05) < Ag max < 0.12(0.08) for an actual intrinsic axial 
ratio of e m i n /e max = 1/3 (2/3). The range for the 68.3% of the 
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Figure 4. The same as in Fig.[5]for e 1 = 1.2, e 2 = 0.8 and e 3 = 1 (T = 0.5). 



as triaxial 




0.2 0.4 0.6 0.8 1 
T 



as prolate 



as oblate 




0.2 0.4 0.6 0.8 
T 



0.2 0.4 0.6 0.8 1 
T 



Figure 5. Inferred axial ratio as a function of the triaxiality parameter T under different hypotheses. The true q max is fixed to 1/3, as shown by the horizontal 
lines. The full line represents the inferred g max averaged over the orientations of the line of sight, assumed to be randomly distributed. The shadowed and 
light shadowed regions contain the 68.3% and the whole range of the inferred values, respectively. q max has been derived under the hypothesis of a triaxial 
ellipsoid aligned along the line of sight, of a prolate and of an oblate ellipsoid in the left, central and right panel, respectively. 




Figure 6. The same as in Fig.HJfor an intrinsic f/max = 2/3. 
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orientations is of 0.01(0.01) < Ag max < 0.23(0.14), whereas the 
total range spans from the true value to a maximum of q max ~ 1, 
i.e an apparently round geometry, for T ~ 0.4 (0.45). 

If we make the hypothesis that the cluster is prolate, the error 
decreases with the cluster effectively approaching a nearly prolate 
shape (T — > 1). For e m i n /e max = 1/3 (2/3), the average devia- 
tion is as large as > —0.09 (—0.11) for T — > 0, and is really small 
(j Ag max | < 0.02) for T > 0.4(0.5). The range for the 68.3% of 
the orientations is of -0.22(-0.22) < Ag max < 0.08(0.03) for 
T = 0, -0.13(-0.08) < Aq max < 0.12(0.09) for T ~ 0.4 (0.6) 
and goes to for T -> 1. The full range is of -0.22(-0.22) < 
Ag max < 0.24(0.15) for T = 0, -0.13(-0.11) < Ag max < 
0.59(0.29) for T ~ 0.4 (0.5) and, as usual, goes to for T -> 1. 

Under the hypothesis of an oblate geometry the situation is 
reversed. For e m i n /e max = 1/3 (2/3), the average deviation is 

small (|Aq max | < 0.02) for T < 0.3 (0.4) and goes to r 0.14 

(—0.13) for T — > 1. The range for the 68.3% of the orienta- 
tions starts from at T = 0, is of -0.09(-0.06) < Ag max < 
0.11(0.08) for T ~ 0.2 (0.2) and goes to -0.22(-0.22) < 
A<? max < -0.03(-0.01) for T = 1. The full range is of 
-0.14(-0.12) < Ag max < 0.59(0.32) for T ~ 0.4 (0.45.) and 
-0.22(-0.22) < Aq max < 0.22(0.3) for T = 1. 



7 BAD IDENTIFICATIONS OF AXIALLY SYMMETRIC 
ELLIPSOIDS 

As seen in Sec. [5] even under the strong assumption that we are 
observing an ellipsoid of revolution, the derivation of the intrinsic 
parameters is in general not unique. In fact, for an axially symmet- 
ric cluster whose size along the line of sight is intermediate with 
respect to its width and the length in the plane of the sky, Eq. {53}, 
we can not establish if it is either prolate or oblate. 



7.1 Prolate as oblate 

A prolate cluster can have the same projected map of an oblate sys- 
tem. An oblate-like deprojection for a prolate ellipsoid is possible 
only when Eq. d52t holds. In terms of the intrinsic parameters of the 
prolate cluster, the inclination angle has to be greater than a given 
threshold, 



cos i < 



VI + 3int ' 



(54) 



if i ^ 7r/4, there is always a possible, but wrong, oblate-like solu- 
tion. As can be seen from Fig. [7] assuming a population of prolate 
ellipsoids randomly oriented, most of them will project in the sky 
in a way that is compatible with an oblate morphology. 

An erroneous hypothesis on the morphology affects the es- 
timate of the intrinsic parameters. Here, let us consider the error 
made when considering as oblate an intrinsically prolate cluster. If 
Eq. J54t holds, the inversion can be performed to infer wrong val- 
ues for qint and i which are still compatible with the observations. 
To find such wrong parameters and infer the error made badly iden- 
tifying a prolate system for an oblate one, we have to substitute in 
Eqs. d50l 15 1 b the expressions of e p and eA for a prolate cluster, 
Eqs. J431 144b . The inferred axial ratio and cosine of the inclination 
angle are 



(55) 



[1- 
and 
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gintl 
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(56) 
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respectively. Results are shown in Fig. [7j In general, the inferred 
axial ratio is less than the actual one. The inferred wrong parame- 
ters are sensitive to the intrinsic axial ratio but depend very weakly 
on the inclination angle. In general the inferred axial ratio increases 
with the true axial ratio. From Eqs. l !54ll55t , we get an upper limit 
for the inferred wrong axial ratio, which is ^ ^/q-mt- For<fc n t < 0.4, 
the inferred value trails the actual one by 0.2; for q ln t > 0.7, the 
gap is reduced to ~ 0.1. Whatever i, an intrinsically pretty elon- 
gated prolate cluster (gi nt < 0.2) can be misidentified for an oblate 
object with a large inclination (cosi < 0.2). On the other hand, a 
nearly round even if slightly prolate object (gi n t > 0.8), can appear 
as an oblate cluster with 0.6 < cos i < 0.7. 



7.2 Oblate as Prolate 

Let us now consider when and how an oblate cluster is projected 
in a way compatible with a prolate geometry. As can be seen from 
Eq. (|47j and Eqs. 0U|49]>, the condition to be fulfilled is 



cosz < 



9int 



1 + qm 



(57) 



Comparing the areas in Fig.|7jand[8]where a double solution is al- 
lowed and assuming a random orientation, we can see how is most 
likely that a prolate objects is interpreted as oblate that the inverse 
case. Under the assumption that an oblate cluster is prolate, one 
will get the wrong intrinsic axial ratio 
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and the wrong inclination (expressed in terms of the cosine) 
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As can be seen from Eqs. f57ll58t . the upper limit for the inferred 
wrong axial ratio is ^qiat, as in the previous case. Results are 
shown in Fig. [8] The inferred axial ratio is less than the actual one, 
i.e. to compensate for the wrong assumption, the cluster is supposed 
to be more elliptical. Errors can be very large. An highly flattened 
oblate object nearly edge-on (Qi nt i; 0.2, cos i < 0.2) can be inter- 
preted as a very elongated prolate cluster nearly face-on. 



8 DISCUSSION 

The combined analysis of X-ray, SZE and GL maps can determine 
the elongation of the cluster along the line of sight, its width and 
length in the plane of the sky and an unbiased estimate of the Hub- 
ble constant. Even if the temperature, metallicity and the density 
profile of the ICM can be determined, the shape and orientation of 
the cluster can not be fully constrained. 

The restrictive assumption of an axially symmetric geometry 
is in general not enough. In fact, a prolate cluster can cast on the 
plane of the sky in the same way of an oblate ellipsoid with different 
inclination and axial ratio and vice versa. Then, in general one must 
assume that the cluster is either prolate or oblate. Even if for some 
particular combinations of orientation and intrinsic ellipticity only 
one axially symmetric solution is possible, there are always some 
more triaxial configurations compatible with data. 
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Figure 7. Intrinsic (wrong) parameters of an intrinsically prolate ellipsoids inferred under the oblate hypothesis as a function of the true intrinsic parameters. 
In the shadow regions, the wrong inversion can not be performed. Left (a) panel: inferred axial ratio g; nt . Thick contour values run from 0.9 to 0.1 with steps 
of 0.1; thin contour values run from 0.09 to 0.01 with steps of 0.01. Right (b) panel: inferred inclination cos i. Contour values run from 0.1 to 0.7 with steps 
of 0.1. 
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Figure 8. Intrinsic (wrong) parameters of an intrinsically oblate ellipsoid inferred under the prolate hypothesis as a function of the true intrinsic parameters. In 
the shadowed regions, the wrong inversion can not be performed. Left (a) panel: inferred axial ratio <fc nt . Contour values run from 0.1 to 1 with steps of 0.1. 
Right (b) panel: inferred inclination cos i. Contour values run from 0.9 to with steps of 0.1. For inferred values smaller than 0.4, contours are pretty close 
to each other. 
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Under the hypothesis that the cluster is prolate or oblate, the 
inversion is likely to estimate the cluster more elliptical than the 
actual value. The main shortcoming is that forcing an axially sym- 
metric geometry to fit an actually triaxial cluster can strongly bias 
the analysis: even the projection of a nearly prolate, but intrinsi- 
cally triaxial cluster can be compatible, for some orientations of 
the system, only with an oblate geometry but not compatible with 
a prolate one and vice-versa. On the other hand, assuming that the 
cluster is triaxial and aligned with the line of sight is a more conser- 
vative approach. The cluster looks rounder and the minor to major 
axis ratio is likely to be overestimated by ~ 0.1, but we are assured 
that the intrinsic geometry is not completely mistaken. 

A very useful additional information for the inversion would 
be the knowledge of the orientation in the plane of the sky of the 
projection of one of the intrinsic axes of the ellipsoid. Whereas 
this can not be obtained by analysing projected maps, other ap- 
proaches could be viable. Unless clusters are axially symmetric 
bodies with internal streaming motions about a fixed symmetry 
axis, line of sight rotational motions along apparent major axis are 
in gener al associated with v elocity gradients along apparent mi- 
nor axes {Contopoulos 1956). This effect can then be studied under 
suitable approximations, i.e. that the velocity field is made up of an 
overall figure rotation about either the shortest or the longe st prin- 
cipal axis together with internal streaming around that axis teinnevl 
Il985h . These theoretical velocity fields can then be projected on to 
the plane of the sky and compared with observations. However, cur- 
rent det ections of velocity gradients in galaxy clusters are still un- 
certain dHwang & Leell2007h ISerencfeOOl and references therein). 

A more direct approach to map out the three dimensional 
structure of a rich cluster is through distances measurements for in- 
dividual members. Distances could b e obtaine d using the method 
of surface brightness fluctuations dMei et all 120071) or empirical 
relationships based on Tully-Fisher or D-a distance indicators 
(Mas ters et afll2006l) . This method is limited to very nearby clus- 
ters, but would give a complete 3-D description of the cluster 
galaxy distribution. 

The theoretical analysis performed in this paper has the advan- 
tage of clearly showing the main degeneracies of the inversion pro- 
cedure. The simple isothermal /3 model was enough to show what 
we can learn on the 3-D structure of galaxy clusters from projected 
maps. Such an analysis avoided some final % 2 fitting procedure to 
some particular simulated map, a method which could miss some 
important facets of the degeneracy question. 

In this paper, I have considered projected X-ray, SZE and GL 
maps. In any case, luminosity and surface brightness observations 
in the optical band would provide the same kind of information. 
Furthermore, if some assumptions are relaxed and we let the galaxy 
distribution have a different intrinsic shape with respect to the ICM, 
we should determine two more intrinsic axial ratios. In any case, 
some other observational constraints could be used in this case. 
Assuming that galaxy and ICM 3-D distributions have the same 
orientation but different ellipticities, then their isophotes w ould be 
misaligned in projection ( Roma nowskv & K ochanek 1998). 

The hydrostatic equilibrium assumption, used in this paper to 
relate GL observations to the other maps, could be substituted by 
other theoretical hypotheses. The assumption of a measurable and 
constant baryon fraction has been suggested as a different way to 
break t he degeneracy between the Hubble constant and the elon- 
gation l lCooravl [l998). However such different theoretical assump- 
tions all aim to break the same degeneracy. If used together they 
would over-constrain some facets of the problem, but would not 



shed light on other undetermined features, such as the orientation 
issue. 
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